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Abstract — In this paper, we present a comprehensive study of 
the monotonicity and log-concavity of the generalized Marcum 
and Nuttall Q functions. More precisely, a simple probabilistic 
method is firstly given to prove the monotonicity of these two 
functions. Then, the log-concavity of the generalized Marcum 
function and its deformations is established with respect to 
each of the three parameters. Since the Nuttall Q— function 
has similar probabilistic interpretations as the generalized Mar- 
cum Q— function, we deduce the log-concavity of the Nuttall 
Q function. By exploiting the log-concavity of these two func- 
tions, we propose new tight lower and upper bounds for the gen- 
eralized Marcum and Nuttall Q functions. Our proposed bounds 
are much tighter than the existing bounds in the literature in most 
of the cases. The relative errors of our proposed bounds converge 
to as b ^ oo. The numerical results show that the absolute 
relative errors of the proposed bounds are less than 5% in most 
of the cases. The proposed bounds can be effectively applied to 
the outage probability analysis of interference-limited systems 
such as cognitive radio and wireless sensor network, in the study 
of error performance of various wireless communication systems 
operating over fading channels and extracting the log-likelihood 
ratio for differential phase-shift keying (DPSK) signals. 

Index Terms — Generalized Marcum Q function, log- 
concavity, monotonicity, Nuttall Q function, tight bounds. 



I. Introduction 

THE generalized Marcum Q— function has a long history 
in the study of target detection probability in radar com- 
munications 1^, (21 ■ Recently, it has gained much attention 
for its important applications in digital communications over 
fading channels, such as the bit error performance analysis 
dealing with differential coherent and non-coherent detection 
1 3 1-| 1 1 1, the efficient performance evaluation of dual selective 
combining diversity over correlated Rayleigh and Nakagami- 
m fading channels (Sj and Il2l - fr3l . the information- theoretic 
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analysis of Rician MIMO channels fT4l-fT5l, the energy 
detection of unknown signals over fading channels 1 16| and the 
primary user detection performance for cognitive radio with 
multiple antennas fVl\. For a > and b > the generalized 
Marcum (J— function of real order > is defined by |18| 

Q,{a,b) = ^j^ ee-'-^I,-i{at)dt, (1) 
Jb 

where is the modified Bessel function of the first kind 
of order v and the right hand side of the above equation is 
replaced by its limiting value when a = 0. If = 1, then 
this reduces to the standard (first-order) Marcum Q— function, 
denoted as Q{a,b). The standard Nuttall (J— function is a 
generalization of the Marcum (J— function, and is defined by 
ESI eq. 86] 

Q^.Aa-b)^ t^'e-*-^h{at)dt, (2) 

Jb 

where b,ix,v > and a > 0. A normalized version of the 
Nuttall Q— function with respect to the parameter a is defined 

as EOI, ED 

Q,,,{a,b) ^ Q^'-^"'^) ^ r ^^e-'^lM)dt. (3) 

In particular, \f ^ = v + \, then Q^,;y(a, 6) reduces to the 
generalized Marcum Q— function of order v + 1, i.e. we 
have Q;y+i,;y(a, 5) = (5y+i(a, 6) for all admissible values of 
a, b and v. The applications for the standard and normalized 
Nuttall (3— function include array processing performance in 
fading channels with co-channel interference |22|, the largest 
eigenvalue of noncentral complex Wishart matrices |23|, per- 
formance analysis of multiple-input-multiple-output (MIMO) 
systems employing multichannel beamforming in arbitrary- 
rank Ricean channels |24| and polynomial approximations 
for evaluation of the average error probability over slow- 
fading channels |25|. Some more applications of the Nuttall 
(3— function can be found in |21 1 and in the references therein. 

However, the precise computations of the generalized Mar- 
cum and Nuttall (J— functions are quite difficult, mainly 
because the modified Bessel function of the first kind 
is involved in the integrands of Q and (|2j. In the last 
few decades, many researchers were working on precise and 
stable numerical calculation algorithms of the generalized 
Marcum and Nuttall (J— functions (see 1261 -lITl | for example). 
While only integer order generalized Marcum and Nuttall 
(3— functions are considered in these papers, appropriate eval- 
uating methods for non-integer order generalized Marcum and 
Nuttall (J— functions are desirable for the performance analysis 
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of some wireless communication systems. Some examples can 
be found in ri2'|-ri3T, [21^ and the references therein. In 
[32|, an exact expression of Q^{a,b) was given for the case 
when v is an odd multiple of 0.5. More compact closed-form 
expressions were proposed in |33| and f2T\. A closed-form 
expression of Q^^^{a^h), where ^^v are odd multiples of 0.5 
and fJ, > I', was also given in 1211 . During the review period 
of this paper, we discovered the papers 1121 and 1131 . in which 
an finite-integral representation of Q^{a,h) with real-order v 
was provided. But this finite-integral representation is still not 
in closed-form. 

To obtain simpler evaluation methods for the generalized 
Marcum and Nuttall Q— functions, one is often willing to 
accept closed-form bounds of the functions if they are tight 
|U4|. In recent years, many new lower and upper bounds were 
proposed as simpler alternative evaluating methods for the 
generahzed Marcum and Nuttall Q-functions f2Tl, f32l-f42l. 
Among them, many tight bounds for the Marcum Q— function 
and the generalized Marcum Q-function were obtained by 
exploiting the bounds of the integrand of ([0 or modifying the 
integral region via a geometric interpretation of the functions 
II34I - I42I . However, these bounds can be only tight for a part 
of the region of h, i.e. either h < a or b > a. 

Recently, the monotonicity of the generalized Marcum and 
Nuttall Q— functions has drawn special interests to provide 
new bounds L21J . [32 1. Different analytical proofs of the 
monotonicity of the function v ^ Q^{a,b) were given in 
1211 . 1431 . However, this result was deduced also with a 
simple probabilistic method in 1441 . where it is also shown 
that the monotonicity of h- > Q^{a,b) is equivalent to the 
monotonicity of the incomplete gamma function ratio, proved 
by Tricomi in 1950 |45|. A summary of the monotonicity 
of Q,y{a,b) with respect to the parameters a and b can 
be found in 1461 p. 451]. The monotonicity of Q^.iy{a.h) on 
H + V, under the requirements that a > 1 and n> v + 1 and 
for constant difference ji — v, was given in 12 II . In this paper, 
a simple probabilistic method is provided which proves in an 
unified way the monotonicity of the generalized Marcum and 
Nuttall Q— functions with respect to different parameters. By 
this method, all the aforementioned results can be obtained, as 
well as a novel result that the normalized and standard Nuttall 
Q— functions are strictly increasing in a. 

The bounds derived via the monotonicity of the generalized 
Marcum and Nuttall functions hold true for all range of b. 
But they are not tight enough in terms of relative errors for 
a part of parameter region. Particularly for the case of large 
b, the relative errors of these bounds are unbounded (for the 
upper bounds) as h approaches infinity. 

In order to get even tighter bounds of the generalized 
Marcum and Nuttall Q— functions, some stronger property of 
these two functions needs to be established. In |4 p. 81], an 
asymptotic formula (but not a bound) for Q,^{a, b) is provided 
when b tends to infinity, given by 

Q,{a,b)^(-\ Q{b-a), (4) 

where Q(-) is the Gaussian Q— function and ^ means 
that these two functions tend to be equal as b increases. 



This impUes that \/Qv-o.5{a, fe)Qiy+o.5(a, b) can be used 
to estimate Qi,{a,b) for integer i' and very large b, and 
the relative error of this estimation converges to when 
b ^ oo. Since several exact expressions of Q^{a,b) for 
odd multiple of 0.5 order v were obtained in 1321 - 1331 
and |21|, ^Qi/-o.5(a, b)Q^+o.5ia, b) can be expressed in 
closed-form. Moreover, our further investigation showed that 
V Qu-o.5{a, b)Qu+o.5{a, b) is also a lower bound of Q^{a, b) 
of integer order v, which is equivalent to 

log Q^(a, 6) > i[logQ^_o.5(a,6) + log(9,.+o.5(a,fe)], 

for integer v. (5) 

This motivated us to start working on the log-concavity of the 
generalized Marcum Q— function and the Nuttall Q— function. 

In 1441 . we proved the strict log-concavity of the functions 
b ^ Qw{cL, b) and b ^ Qu{a, Vb) for all > 1 and a > 0. 
Moreover, based on some preliminary results, we conjectured 
that the function ly i— > Q^{a,b) is strictly log-concave on 
(0, oo) for all a > and & > 0. In this paper, we are 
able to verify the above conjecture on [l,oo). We present a 
comprehensive study on the log-concavity of the generalized 
Marcum and Nuttall Q— functions. New tight bounds for 
real order generalized Marcum and Nuttall Q— functions are 
proposed based on the log-concavity of these functions. Our 
proposed bounds are much tighter than the existing bounds 
in the literature in most of the cases. They involve only 
exponential function and the erfc function, and therefore 
can be computed very efficiently. The relative errors of the 
proposed bounds converge to when b ^ oo. The numerical 
results show that the absolute relative errors of the proposed 
bounds are less than 5% in most of the cases. To the extent 
of the authors' knowledge, our bounds for the generalized 
Marcum and Nuttall Q— functions are the first bounds with 
such tightness in terms of relative errors on the whole region 

of & e (o,oo). 

The detailed content is as follows: 

In Sectionln] we first review the probabilistic interpretations 
of the generalized Marcum and Nuttall Q— functions related 
to the non-central chi and chi-square distribution, which form 
the basis of this paper 

Then, in Section |llll we present a simple probabilistic 
method which proves in an unified way the monotonicity of 
the generalized Marcum and Nuttall Q— functions with respect 
to different parameters. 

In Section lTvl we first recall a mathematical concept named 
total positivity, which plays an important role in the following 
proofs. In complement to the results of Theorem 2.7 in |44|, 
we prove that the functions 6 i-^ 1 — Qv{a, Vb), b t-^ 1 — 
Qv{a, b) are log-concave on (0, oo) for all > 1 and v >2>/2, 
respectively. Then, we show that the functions v ^ Q^{a,b) 
and ly 1-^ 1 — Qv{a,b) are log-concave on [1,00), while the 
functions v ^ Qi/(0, b), v 1— > Qu{a, b) — Q„{0, b) and 1/ 1-^ 1 — 
(5i,(0, 6) are log-concave on (0, 00), for all admissible values 
of a and b. We also prove that the functions a 1-^ (v^i b), 
a ^ \ — Qu (-v/flj b) and a ^ \ — {a, b) are log-concave 
on (0, 00). Some remarks and conjectures are also provided. 

In view of the close relationship between the Nuttall 
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Q— function and the generalized Marcum Q— function, the log- 
concavity of the Nuttall Q— function and its deformations are 
established in Section IVl 

Closed-form lower and upper bounds are proposed for the 
generalized Marcum and Nuttall Q— functions in Section IVII 
Numerical results and rigorous analysis are also provided to 
justify the tightness of the bounds. 

Some applications of the theoretical results are provided in 
Section IVTIl 

Finally, the conclusion of this paper is given in Section rvilll 
For ease of latter use, we define the following notations: 
The probability density function and characteristic function 
of a random variable X are denoted hy x fx (x) and 1 1-^ 
(px{t), respectively. We use Ex{X) to denote the expectation 
of X, and Ey{Y\X — x) to denote the conditional expectation 
of Y given at the value X = x. The indicator function of a 
set S is defined as 



1, if a: e 5, 
0, otherwise. 



(6) 



We use M to denote the set of real numbers, N to denote the 
set of positive integers. 

When we discuss the log-concavity of a function /, it is 
convenient to allow / to take on the value zero, in which case 
we take log /(a;) = —oo |47l, f48l. By this means, we have 
that the indicator function Cs of a convex set S is log-concave. 

II. Probabilistic interpretations of the 

GENERALIZED MARCUM AND NUTTALL Q-FUNCTIONS 

The generalized Marcum Q— function Qi,{a,b) has an im- 
portant interpretation in probability theory: it is the com- 
plementary cumulative distribution function (ccdf) or relia- 
bility function of the normalized non-central chi-square dis- 
tribution with degrees of freedom |4] p. 82]. For this, 
let Xi,X2,-.-,X^ be random variables that are normally 
distributed with unit variance and nonzero mean 7,, where 
i £ {1,2, . . . , I'}. It is known that the random variable 
Xf + X2 + • • • + X^ has the non-central chi-square distribu- 
tion with I' degrees of freedom and non-centrality parameter 
+ 7^. The probability density function 



7l 



A = jt 
(pdf) : (0,(X3) (0,00 

distribution is defined as 



of the non-central chi-square 



42 (x) = 2-''/2e-(^+^)/2 



fe>0 



i//2+fc-l 



(A/4)* 



r{u/2 + k)k\ 



-{x+\)/2 



(l^j (VA^) , (7) 

where F is the Euler gamma function [W eq. (6.1.1)]. Al- 
though V is an integer in our description of the non-central 
chi-square distribution above, it is known that the distribution, 
defined by 0, is a proper distribution for any positive v 
II46I p. 436] and |44|. Consequently, the generalized Marcum 
Q— function can be expressed as 



2 \aJ 



i>/2-l/2 



.-{x+a)/2j^_^ (V^) dx. (8) 



In view of the similarity of the definitions of the generalized 
Marcum and Nuttall Q— functions ([0 and Q, one can find 
that the normalized Nuttall Q— function is actually the upper- 
side partial (truncated) moment of the normalized non-central 
chi-square distribution, which is given by 



{x)dx 



00 ^(f,-l)/2 

2a''/2 " 



(v^) dx. 



(9) 



If A = 0, the above distribution reduces to the classical 
(central) chi-square distribution, whose pdf is 



X' 



iy/2-1 ~x/2 



'x,A I Jx,..„y I 2''/2r(i//2) 
Consequently, (|8j reduces to 



(10) 



(0, 



f^2 (x)dx = / — — dx. (11) 



We note that if the random variables , and - , 
which have non-central chi-square distribution with i/i, 1/2 
degrees of freedom, respectively, and non-centrality param- 
eters Ai, A2, respectively, are independent, then the random 
variable xt^ Xi + X2 ^ non-central chi-square distri- 
bution with vi + V2 degrees of freedom and non-centrality 
parameter Ai + A2. Therefore, the non-central distributions 
is reproductive under convolution [46J. This can be explained 
by the characteristic functions of the non-central and central 
chi-square distributions, which are given by 



and 



^^2 f^)=eT^(l-2it)-''/2, 



(12) 



(13) 



It may be a little surprising that the degrees of freedom v 
can be in J12> . and the classical non-central distribution 
reduces to non-central x^ distributions of zero degrees of 
freedom, i.e. Xq ^ |50|. Xo a '^^'^ ^1^° defined as a mixture 
of the distributions 0, x^, xli ■ ■ ■ with Poisson weights. 
The above mentioned reproductive property holds true for 
non-central distributions of zero degrees of freedom, 
i-e- xlx, + Xlx2 ~ Xo,Ai+A2' xlx + xl - xlx and 
Xo,Ai +X^,A2 ~ X^.Ai+A2- However, the random variable Xo,a 
does not possess a probability density function, because of a 
discrete mass of e^^^ at the zero point, and the positive part 
of this distribution has a density function /^2 ^ , by excluding 
the probability at a; = 0, which can be expressed as 



1 

—I 

X 

I /\ \ 1/2 

2 



fc!(fc- 1) 



-(^+^)/2/lf^/A^' 



(14) 
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The ccdf of Xq \ is 

P {xlx > b) 

\ l\{\)''\-^-^h(y^dx^ if 6>o, ^^^^ 

1, if & < 0, 

which has a discontinuity point at 6 = 0. 

There are also alternative probabilistic interpretations of the 
generalized Marcum and Nuttall Q— functions related to the 
non-central chi distribution. It is known that if Yi, I2, ■ ■ • , i^i/ 
are random variables that are normally distributed with unit 
variance and nonzero mean /iti, /i2, . . . , /i,/, then the random 
variable ^Y^ + . . . + Y^]^^^ has the non-central chi distiT- 
bution with ly degrees of freedom and non-centrality parameter 
T= [/i?+A'i + --- + M^]'/'-Thepdf/x_ : (0,00)^ (0,00) 
of the non-central chi distribution |46| is defined as 



k>0 



2k 



r{iy/2 + k)ki 



re 



(16) 



where v can be also an arbitrary positive real number There- 
fore, the generalized Marcum Q— function and normalized 
Nuttall Q— function can be expressed as 



lb 



(17) 



and 



poo 

Q^,,{a,b)^ x^-''-^fx,^^^,^Sx)dx, (18) 

Jb 

respectively. If r = 0, the above distribution reduces to the 
classical chi distribution with pdf fx„a ■ (Oi 00) ~* (0, 00), 
defined by 



Therefore, ilH reduces to 



Qu (0,6) 



b Jb 



00 ^2u-l^-x^/2 



(19) 



dx. (20) 



Some more probabilistic interpretations of (a, 5) can be 
found in C0|, |32| and l5Tl-l52l. 

The chi and chi-square distributions have close connection 
to a number of useful distributions in digital communications, 
such as Rayleigh distribution, Rician distribution, Nakagami- 
m distribution, generalized Rayleigh and Rician distribution, 
bivariate Rayleigh and Nakagami-m distribution, exponential 
distribution, gamma distribution, etc. OJ-H, |53|-|54|. This is 
exactly the reason for the vast applications of the generalized 
Marcum and Nuttall Q— functions in wireless communications. 

'There is a mistake in the original foimula given in Siegel's paper 1501 . 
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Fig. 1. A geometrical intuition of our proof idea of Lemma ITI 



III. The monotonicity of the generalized Marcum 
AND Nuttall Q-functions 

The monotonicity of the generalized Marcum and Nuttall 
functions has been discussed via a number of analytical 
proofs 1211 . 1321 . 1431 . 1461 . However, these proofs did not 
provide much insight about the functions, and each of proof 
could only show the monotonicity with respect to one param- 
eter. 

In this section, we utilize the powerful probabilistic inter- 
pretations of the generalized Marcum and Nuttall functions 
that we have discussed earher, and present a simple probabilis- 
tic method to solve in an unified way the monotonicity of these 
functions with respect to different parameters. A geometrical 
intuition of our method is also provided to facilitate the 
understanding of the proof. 

In part (a) of Theorem 3.1 in |44|, we have used a simple 
probabilistic result for independent random variables X and 
Y with non-negative support, i.e. 



P{X + Y >b)> P{X > 6), for 5 > 0, 



(21) 



to prove that i/ 1-^ Qi,{a,b) is strictly increasing for a > 
and 6 > 0. 

The inequality J21I I can be explained simply in a geometrical 
way, as illustrated in Fig. [J We consider the first quadrant 
in the coordinate plane of x and y, where A is the region 
{{x,y)\x > b,y > 0} and B is the region {{x,y)\x + y > 
b,0 < X < b,y > 0}. Then, P{X + F > 6) is the sum of 
the probabilistic integrations on the regions A and B, while 
P{X > b) is the probabilistic integration only on the region 
A. By this, J2U is expected. 

In order to solve in an unified way the monotonicity of 
the generalized Marcum and Nuttall Q— functions, we now 
generalize the inequality i21\ by inducing an non-decreasing 
positive weight function in the probabilistic integrations. Our 
new result is stated as 

Lemma 1. Let X and Y be independent random variables 
with non-negative support and let g : (0, 00) — > (0, 00) 
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be an non-decreasing positive function. Further, let fx and 
fx+Y be the pdfs of the random variables X and X + Y, 
respectively. Let Fy be the cumulative distribution function 
( cdf) of the random variable Y, which may not has a pdf if 
Fy is discontinuous. Then, the following inequality is true for 
each b > 

f g{t)fx+y{t)dt> f g{t)fx{t)dt, (22) 

Jb Jb 

if Fy{0) < 1 and the integrals exist. 

The proof of Lemma [2 is given in Appendix IaI which also 
fits the geometrical intuition of Fig. ^ 

With this Lemma, we now can prove the strict monotonicity 
of the generalized Marcum and Nuttall Q— functions. 

We first consider the case when g{x) = 1 for all x e (0, oo). 
If X ^ xLi,a and Y ~ xL,, then clearly X + Y ^ 
^2(1/1+1/2) a- Hence, by Lemmafjwe easily have that (see also 
Theorem 3.1 in [44]) ly Q,^{^/a,^/b) is strictly increasing 
on (0, 00) for each a > and & > 0, i.e. we have 

(^y/a, Vbj > Q^^ (^y/a, Vbj (23) 

for all 1^1, z^2, > and a > 0. On the other hand if X ^ 

X2v,ai aid Y ~ Xo,a2 ' ^^^^ ^° ^^^^ X + Y ^ 

xiiyai+a,- ^c note that Xo a, possesses a discontinuous cdf, 
and therefore has no pdf. Consequently, in view of Lemma 
[2 we obtain that the function a ^ Qu{\fa,"s/h) is strictly 
increasing on [0, 00) for each z^, 6 > 0, i.e. 

Qu (Vai +a2, VPj > Qu (Vai, ^) (24) 

for all fli > and 02, i^, 6 > 0. 

Finally, let g{x) = x^/^ for all x € {0,oo) and ^ > 0. A 
similar argument as we presented above yield the following 
inequalities for the normalized Nuttall Q— function 

(25) 

and 

Q^L+u.v-l {\/ai + 0.2, Vtj > Qt_i+iy,u-l (^y/oi, \/Pj (26) 

where I'l, > 1, /i > and ai, a2, j^2, a, & > 0. Moreover, 
it is easy to verify that the generalized Marcum and Nuttall 
Q— functions are strictly decreasing with respect to h on 
[0, 00). These results are summarized in the following theorem. 

Theorem 1. The following assertions are true: 

(a) The generalized Marcum Q— function Qv{a,b) is strictly 
increasing in v and a for all a > and b,u > 0, and is 
strictly decreasing in b for all a,b > and v > 0. 

(b) The normalized Nuttall Q— function Q^+i,_(_i ,y(a, 6) is 
strictly increasing in v and a for all jJ. > Q and 
a,b > 0, and is strictly decreasing in b for all J^, /i > 
and a > 0. 

(c) The standard Nuttall Q— function (5^_|_iy_|_i^iy(a, 5) is 
strictly increasing in v for all a > 1, v, fi > and 6 > 0, 
is strictly increasing in a for all v, fJ. > and a, 6 > 0, 
and is strictly decreasing in b for all b,iy, ^ > and 
a > 0. 



-1.5 



-1.6 - 




-2.3 1 • • • • • ' • 1 

1 1.5 2 2.5 3 3.5 4 4.5 5 

v=x 

Fig. 2. The logaiithm of the pdf of non-centi'al chi-square distribution, i.e. 
log (x), when A = 2 and u = x. 



-1.4 



-1.6 



-1.8 



-2 



-2.2 



-2.4 



-2.6 




Fig. 3. The logarithm of the pdf of non-central chi-square distribution, i.e. 
log (x), when v = 2 and X = x. 

Therefore, we have obtained the monotonicity of the gener- 
alized Marcum and Nuttall Q— functions with respect to vari- 
ables a, b and with an unified and rather simple probabilistic 
method. To the best of the authors' knowledge, the mono- 
tonicity of the standard and normalized Nuttall functions 
with respect to the non-centrality parameter a has not been 
discussed in the literature. 

IV. The log-concavity of the generalized Marcum 

Q-FUNCTION 

In II44I . Sun and Baricz conjectured that the function f 1-^ 
Q^{a,b) is strictly log-concave on (0,00) for all a > and 
6 > 0. In order to prove this, we have tried several methods. 
One idea is to consider the generalized Marcum Q— function 
as a function of two variables and to prove that it is log- 
concave as a function of these two variables. It is known that 
the pdf of the non-central chi-square distribution is log- 

concave with respect to h' and is also log-concave with respect 
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to a under some requirements II55I . However, these results do 
not help too much in the study of the log-concavity of the 
functions v i-^ Q^{a,b) and a ^ Q^{\/a, b), because we may 
need the function {v, x) i— > ^ (x) or (A, x) i-^ f^2 ^ (x) 
to be log-concave on two variables |47 p. 106]', which is 
not true. In Fig. |2] and |3] we show the numerical results of 
the functions x i-^ log ^ (x) and x i-^ log (a;), from 
which one can understand that (i^, x) i— > f^2 (x) or (A, x) i-^ 
f^2 ^ (x) can not be log-concave. 

Fortunately, we have found a powerful mathematical con- 
cept, named total positivity, which can help us to conquer 
this difficulty. For the reader's convenience, we first offer a 
brief introduction to total positivity, which is necessary for 
our proofs. More details about total positivity can be found in 
Karlin's monograph 1571 . 

A. A brief introduction to total positivity 

The definition of totally positive function fST* p. 11] is the 
following: 

Definition 1. A function f{x,y) of two real variables, x 
ranging over Ai and y ranging over A2, is said to be 
totally positive of order r (TPr), if for all Xi G Ai and 
yi € A2, i G {1, 2, . . . , m} such that xi < X2 < ■ ■ ■ < Xm, 
2/1 < 2/2 < • • • < 2/mj where m G {1, 2, . . . , r}, we have the 
inequalities 



f 



•^1 ) •^2 ) ■ • ■ 7 
2/1,2/2, • ■ • ,2/rr 

f{xi,yi) f{xi,y2) 



f{xi,ym) 



> 0. (27) 



f{xm,yi) /(a;™, 2/2) •■■ f{xm,ym 

Either of Ai and A2 can be an interval of real line, or a 
sequence of countable discrete points. 

A related concept to total positivity is the sign reverse 
regularity. 

Definition 2. A function f{x,y) is said to be sign reverse 
regular of order k (SRRk), if for every xi < X2 < ■ ■ ■ < Xm 
and 2/1 < 2/2 < ■ ■ ■ < 2/m, where m G {1, 2, . . . , k}, the sign 
of 



f 



2/1,2/2, 



,2/n 



(28) 



(_l)™(™-l)/2 



If a function f{x, y) is said to be TP^o, if it is TPk for all 
fc G N. Similarly, it is said to be SRRoo, if it is SRRk for all 
fc G N. 

Totally positive functions can be generated by the following 
composition law, which expresses a continuous analogue for 
evaluating determinants arising in matrix multiplication fSTl 
pp. 16-17]. 

'if f(x, z) is log-convex in z for eacli x £ C, then g(z) = f(x, z)dx 
is log-convex, where C is an arbitrary set of x 147.1 p. 106] and |56 p. 9]. 
However, a similar statement of this does not hold true for log-concavity. 
Instead, the following statement is true for log-concavity: If /(x, 2) : M" X 
R™ 1-^ R is log-concave in x an d z, then g( z) = J fix, z)dx is a log- 
concave function of 2 on R™ El P- 106] and ^ Theorem 3]. 



Lemma 2. Ifr{x,w) — Jp{x,t)q{t,w)da{t) and the integral 
converges absolutely, then 



Xi,X2, ■ 
Wi,W2, ■ 

P 



■■,Xk 
■■,Wk 

Xi,X2, ■ 
tl, t2, ■ 



tl<t2<---<tk 

tl,t2, ■ ■ ■ jtk 
Wi,W2,...,Wk 



d(j{ti)dcf{t2)...da{tk). (29) 



Here, dij{t) denotes a sigma-finite measure defined on t G A. 
When A consists of a discrete set, the integral is interpreted 
as a sum. 

The above result has an immediate corollary, expressed as 
follows J57| 



Lemma 3. If p{x, t) is TPm and q{t, w) is TPn, then 
r{x,w) = J p{x,t)q{t,w)da{t) is TP^i^^j^ nj, provided that 
a(t) is a regular finite measure; when p{x,t) is SRRm and 
q{t,w) is TPn, then r{x,w) is SRRaiin{m.n}- 

Analogously, as in Lemma |2l when t is chosen from a 
discrete set, the above integral is interpreted as a sum. 

The following lemma states another important composition 
law of total positive functions 1571 p. 130]. 

Lemma 4. Let f{x + y) be SRRr for x,y > 0. Suppose that 
(j){t, x) is TPr for t,x > 0, and satisfies 



(t){t + s, x) 



If c is defined by 



c{t) 



cj){t,Ocl){s,x-Odt (30) 



(31) 



then c{t + s) is SRRr for t, s > 0. 

We note that formula ( I30l l expresses exactly the reproductive 
property of (j>{t, x) under convolution. 

Some other results which are useful in the sequel are 
included in the following lemma. The proof of parts (a)-(d) 
can be found in 1571 and part (e) is stated in 147. p. 79]. 

Lemma 5. The following assertions are true: 

(a) f{x, y) = xy is TP2 for Ai = (0, 00), A2 = (-00, 00). 

(b) If f{x,y) is TPr, let (j){x),ip{y) maintain the same 
constant sign on Ai and A2, respectively, then 
cj){x)ip{y)J{x,y) is also TPr. 

(c) Let f be a strictly positive second order differentiable 
function. Then {x,y) t-^ f{x — y) is TP2, if and only if 
X I— !■ f{x) is log-concave on its domain. 

(d) {x,y) I— > f{x + y) is SRR2, if and only if x ^ f{x) is 
also log-concave on its domain. 

(e) Let / : M" M, A G M"^", and b G R", and consider 
the function g : M™ i— )■ M, defined by 



g{x) = f{Ax + b), 



(32) 



with domg — {x \ Ax + 6 G domf}. If f is log-concave, 
then g is log-concave too. 
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B. The log-concavity of Q^{a,b) with respect to b 

Recently, we proved in 1441 that the functions b i-^ 
Q^{a,Vb) and 6 i-^ Qiy{a,b) are both strictly log-concave 
on (0, oo) for all a > and ly > 1. Now we give the 
complementary results as follows. 

Theorem 2. If a> 0, then the functions 6 i— > 1 — Qy{a, ^/b) 
and 6 I— > 1 — Q^{a, b) are log-concave on (0, oo) for all v > 1 
and V > 3/2, respectively. 

Proof: The function 1 — Q^{a,\/b) is the cdf of non- 
central chi-square distribution, given by 

l-g,(a,V6)= / 42 {x)dx. (33) 
Jo 

We know that the function x ^ f^2 (x) is log-concave on 
(0, 00) for all a > and v > 1 \55i . On the other hand it is 
known that if the probabihty density function is log-concave, 
the cumulative distribution function is also log-concave (see 
i58J and the references therein). Therefore, b ^ 1 — Qy{a, Vb) 
is log-concave too on (0, 00) for all a > and v >1. 

We have proved in [441 that the pdf of non-central chi 
distribution b i~> /^^^ ^ (6) is log-concave on (0, 00) for all 
a > and v > 3/2. When a = 0, the pdf of central chi 
distribution b 1-^ /x2i-(^)i defined by jl9> , is log-concave on 
(0, 00) for all > 1/2. From ^ and (|20|i, we have 

1-Q,(a,&) = / fx2.A^)dx. (34) 
Jo 

By using the same method as above, we conclude that the 
function 6 i-^ 1 — Qu{a,b) is log-concave on (0, 00) for all 
a > and v > 3/2, and the proof is complete. ■ 

We have shown the log-concavity of the Marcum Q— 
function and its deformations with respect to b for not quite 
small t/. Now we discuss the cases of small which is 
suggested by our numerical results. 

Remark 1: It was shown that b i— > (5^(0, y/b) is log-concave 
for V e [1, 00) and log-convex for v G (0, 1] in |59|. On the 
other hand, our numerical results show that b Qu{a, Vb) is 
neither log-convex nor log-concave if u E (0, 1] and a > 0. 
Therefore, our previous result that b t-^ Qu{a, Vb) is strictly 
log-concave on (0,00) for all a > and 1/ > 1 is sharp. 

On the other hand, it was shown that the pdf of chi 
distribution is log-concave for ly > 1, which indicates that 
b 1-^ (5j/(0,6) is log-concave on (0,oo) for all ly > 1/2 1591 . 
our numerical results suggest that b 1-^ Q^, (a, b) is log-concave 
on (0, 00) for all a > and i/ > 1/2. Therefore, we conjecture 
that 

Conjecture 1. The function b 1— > Q^, (a, 5) is log-concave on 
(0, 00) for all a > and v e [1/2, 1]. 

A special case of Conjecture^is when v = 1/2. In this case, 
the conjecture reduces to 6 1-^ {a,b) = Q{b-\-a) + Q{b—a) 
is log-concave on (0, 00) for a > 0. 

Remark 2: In |59|, Bagnoli and Bergstrom proved that the 
cdf of gamma distribution is a log-concave function. Since the 
chi-square distribution with ly degrees of freedom is a gamma 
distribution with parameter j//2, they derived that 6 1 — 



Qu{0, Vb) is log-concave for all > 1591 p. 16]. Moreover, 
using the same proof idea for the case of chi distribution, one 
can show that b 1-^ 1 — Qiy(0, b) is log-concave for all ly > 0. 
We now provide a conjecture for the case a > 0, stated as 

Conjecture 2. If a > 0, then the functions b i~> 1 — Q^(a, Vb) 
and 6 I— > 1 — Q^{a, b) are log-concave on (0, 00) for all v G 
[0,1) and V € [1/2,3/2), respectively. 

We note that our numerical results suggest that b ^ \ — 
Qq ^a, Vb^ is log-concave on (0, 00) while b ^ 1 — (a, b) 
is not log-concave for v E (0, 1/2). 

C. The log-concavity of Qi^{a,b) with respect to v 

In 1 44 1, we deduced a Turan type inequality, which is 
interesting in its own right 1551 . given as 

Ql+i{a, b) > Q,{a, b)Q,+2{a, b). (35) 

It is known that the integrand of QviV^-i Vb) as a function 
of V, i.e. V ^ (x) is log-concave on (0,cxd) 1551 . 

Thus, the above Turan type inequality suggests that this log- 
concavity property remains true after integration, of course 
with some assumptions on parameters. Taking into account 
this observation in |44| we conjectured that the function 
ly Qi^{a, b) is actually strictly log-concave on (0, 00) for all 
a > and b > 0. Now, we are able to verify this conjecture 
for v >1. The case v E (0, 1] remains open. 

Theorem 3. The following assertions are true: 

(a) The function v 1— > Q^, (0, b) is log-concave on (0, cxd) for 
all 5 > 0. 

(b) The function v i— )■ Qyia^b) is log-concave on [l,oo)/or 
all a,b > 0. 

(c) The function v 1-^ Q,y{a, b) — Qi,(0, b) is log-concave on 
(0, c») for all a.b > 0. 

(d) The function ly t-^ 1 — Q,y{0, b) is log-concave on [0, 00) 
for all b> 0. 

(e) The function v i-^ 1 — Q^{a, b) is log-concave on [1, 00) 
for all a > and b > 0. 

The proof of Theorem |3] is provided in Appendix 151 
Remark 3: It is worth to mention here that very recently 
Alzer and Baricz |60| by using an interesting idea of Alzer |6^ 
proved that for all 6 > fixed the function v 1-^ Qu{0, V^b) is 
log-concave on (0,oo). Clearly this implies that for all & > 
the function ly Qi/(0, b) is log-concave too on (0, 00), which 
is exactly the statement of part (a) of Theorem|3l However, in 
Appendix |b1 we present a completely different proof for this 
part. 

Moreover, we note that Merkle |62| based on a Turan- 
type inequality involving the incomplete gamma function 
conjectured that the function v ^ 1 — Qy(0, V26) is log- 
concave on (0, 00) for all 6 > 0. A proof of this conjecture can 
be found in Alzer's paper 161]. Clearly, part (d) of Theorem|3l 
verifies also this conjecture, and the proof given in Appendix 
lEIis completely different than in |61|. 

Remark 4: Parts (b) and (e) of Theorem|3]do not include the 
region v E (0, 1]. Therefore, we have the following conjecture 
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Conjecture 3. If a > and b > 0, then the functions v s- 
Quia, b) and v ^-^ 1 ~ Q^{a, b) are log-concave on (0, 1]. 

It is convenient to understand why our proof fails in these 
cases. The key proof idea of Theorem |3l is to use Lemma 
|4l which requires the concerned function with a lower order 
V to be log-concave in b. We mentioned in Remark 1 that 
our numerical results show that b ^ Quia, V6) is neither log- 
concave nor log-convex for ly e (0, 1]. Therefore, our proposed 
method can not be used to prove that v i— > Qy{a,b) is log- 
concave on (0, 1]. Some other methods are needed to prove 
our conjecture. 

On the other hand, we have mentioned in Remark 2 that 
our numerical results suggest b ^ 1 — Qq ^a, \/b^ to be log- 
concave on (0,oo). If one can prove this, our proof method 
can derive that Part (e) of Theorem |3] is true on (0, c»). 

D. The log-concavity of Qu{a, b) with respect to a 

Now, we study the log-concavity of the generalized Marcum 
Q— function with respect to a. It is interesting that the log- 
concavity holds true for the following integral of (x) on 
an interval {c,d) C (0,oo). 

d 

Lemma 6. The function a J f^2 {x)dx is log-concave 

c 

on [0, oo) for all v > and (c, d) C (0, oo). 

The proof of Lemma |6l is given in Appendix [C] 
Our main result of this subsection, which is an immediate 
application of Lemma |6l reads as follows. 

Theorem 4. Let v > 0. Then the following assertions are 
true: 

(a) The function a i— > Qiyiy^jb) is log-concave on [0, oo) 
for all b>0. 

(b) The function a 1 — (5j/(y^, b) is log-concave on [0, oo) 
for all 6 > 0. 

(c) The function a 1 — Qu{a, b) is log-concave on [0, oo) 
for all b> 0. 

Proof: (a) & (b) Using Lemma |6l for the case when c 
tends to b^, and d tends to oo we obtain the result of part (a). 
Similarly, substituting c with and d with b^ in Lemma |6l 
part (b) is also proved. 

(c) It is known that if a positive function / is log-concave 
and decreasing, and g is convex, then the composite function 
f o g is log-concave too 1471 p. 84]. We choose /(a) = 1 — 
Qv{-\/a,b) and g{a) = a?. We known that a y-* f{a) is log- 
concave and strictly decreasing, a i-^ g{a) is convex, thus 
a 1-^ (/ o g){a) = 1 — Qu{a,b) is also log-concave, which 
completes the proof of pait (c). ■ 

V. The log-concavity of the standard and 

NORMALIZED NUTTALL Q-FUNCTIONS 

The standard and normalized Nuttall Q— functions have sim- 
ilar probabilistic interpretations with the generalized Marcum 
Q— function. Therefore, we can establish the log-concavity of 
the Nuttall Q— function and its deformations similarly as in 
the previous section. 



g 

A. The log-concavity of the Nuttall Q— function with respect 
to b 

The following theorem presents the log-concavity of the 
Nuttall Q— function with respect to b, which is similar with 
the results of i44i Theorem 2.7] and Theorem|2lof the present 
paper 

Theorem 5. Let 6, /i, > and a > 0. Then the following as- 
sertions about the normalized Nuttall Q-function, Qfi^u {a, b), 
are true: 

(a) The function b t-^ Qfj,.,^{a, Vb) is log-concave on [0,oo) 
for all fj, + 1/ > 1. 

(b) The function b i—s- Q^.j/ (a, b) is log-concave on [0, oo) 
for all fj, + v > 1. 

(c) The function b -Ex2(„+i) ^i^'^^"^^) ^ Q^i,u{a,Vb) is 
log-concave on (0, oo) for all ii + i> > 1. 

(d) The function b i-> Ex^^^_^_^^^^{X^'-•'-'^) - Q^,^u{a,b) is 
log-concave on (0,oo) /or all fi > 1, v > 1/2. 

The proof of Theorem |5] is given in Appendix |D] These 
results can be rewritten in terms of the standard Nuttall 
Q— function as follows. 

Corollary 1. Let b, > and a > 0. Then the following 
assertions about the standard Nuttall Q-function, Qfi^v ifl, b), 
are true: 

(a) The function b i— > Qf^^uia, Vb) is log-concave on [0,oo) 
for all -\- V > 1. 

(b) The function b Qii,iy (a, b) is log-concave on [0, oo) 
for all fj, + v > 1. 

(c) The function b ^ a^^^xsi^+D.a (^^"''"^) - Qt^-^ia, Vb) 
is log-concave on (0, oo) for all fi + > 1. 

(d) The function b i-^ a'^^x2(,.+i) o {X'^^"^^) — Qfi,v (a, b) is 
log-concave on (0, oo) for all fi > 1, i> > 1/2. 

B. The log-concavity of the Nuttall Q— function with respect 
to the order fi and v 

In this subsection, we consider the log-concavity of the 
Nuttall Q— function with respect to ^ for two different cases 
of the order ^ and v: (1) ^ — v \s fixed, (2) v is fixed. The 
derived results for these two cases are quite different. 

We first study the log-concavity of the Nuttall 
Q— function when ji — v \s fixed. Consider the integral 
of a;(^"'^"i'/2/^2 ^ (x) on an interval (c, d) C (0, oo), the 
following result is true. 

d 

Lemma 7. The function i/ 1— > f x^'^^"^^)/^ f^2 {x)dx is 

C 

log-concave on [0,oo) for fixed /i — > 1, (c, rf) C (0,oo) 
and a > 0. 

The proof of Lemma is given in Appendix |E| By an 
immediate application of Lemma we obtain the following 
result. 

Theorem 6. Let a > and — v > 1 be fixed. Then the 
following assertions about the normalized Nuttall Q— function, 
Qfj.,!^ (flj b), are true: 

(a) The function v ^ Qt^.u (o, b) is log-concave on [0, oo) 
for all b>0. 
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(b) The function v ^ ^x^u^i) ^ ^) ~ Qi-i,}^{a,b) is 

log-concave on [0, oo) for all b > 0. 

Proof: Substituting a with a^, c with and d with 
oo in Lemma we obtain the result of part (a). Similarly, 
substituting a with a^, c with 0, and d with 6^ in Lemma 
part (b) is also proved. ■ 
Since a'' is log-linear in ly, we obtain a similar result for 
the standard Nuttall Q— function. 

Corollary 2. Let a > and fJ- — ly > I be fixed. Then the 
following assertions about the standard Nuttall Q— function, 
Qii.v (i, b), are true: 

(a) The function v i—s- Q^^^v (a, fe) is log-concave on [0, cxd) 
for all & > 0. 

(b) The function v ^ o-''E^2(^+i) a {X'^^"^^) — Q^,i, {a, b) is 
log-concave on [0, cxd) for all b > 0. 

When > is fixed, we may expect the Nuttall Q— function 
to be log-concave in fi. However, the Nuttall Q— function is 
actually log-convex in /i in this case. 

Theorem 7. Let v > be fixed and a > 0. Then the fol- 
lowing assertions about the normalized and standard Nuttall 
Q— function, i.e. Q^^u (a, 6) and Q^.i/ (a, 6), are true: 

(a) The functions /i > Q^,v {o,,b) and /i i-^ Qp,,u (o-ib) are 
log-convex on [0, oo) for b > 0. 

(b) The functions /i i-^ ^x2(^+i) „ '^^^)"~ i'^^b) and 
II |—^ a''-E^2(..+i).a {X'^^"^ ) ~ Q^^i, (a, b) are log-convex 
on [0, oo) for 6 > 0. 

Proof (a) We know that the integrand of (|3} is log-linear 
in ji, and hence is log-convex in for v.^i > {) and a > 0. 
Moreover, if function /(x, y) is log-convex in x for each y e 
C, then 

(36) 



g{x) = / f{x,y)dy 
Jc 

is also log-convex i47l p. 106]. Therefore, Q^,^ (a, 6) 

is log-convex on [0,oo). By multiplying a'^ , we obtain that 

^ Qn,i' (flj is also log-convex on [0, oo) and part (a) is 
proved. 

(b) It is known that 



^X2(.+ l),a(^ 



Q^.^(a,6)= / -^e ' 2" I^{at)dt, 
Jo 

(37) 

which has the same integrand with (|3j, therefore the results 
of part (b) can be proved with the same argument. ■ 



The proof of Lemma|8]is given in Appendix|F| It is actually 
quite similar with that of Lemma |6l Substituting the integral 
limits in Lemma [Sj we obtain the following result. 

Theorem 8. Let fijV > and — v > 1 be fixed. Then the 
following assertions about the normalized Nuttall Q-function, 

Qfj,^u (a, b), are true: 

(a) The function a Q^.u {y/o-, b) is log-concave on (0, 00) 
for b>0. 

(h) The function a ^ E^2 ^ {X^^'-•'-^^/^)-Qf,^^ {^,b) 
is log-concave on (0, 00) for b > 0. 

Proof: (a) & (b) Tending with c to b^, and with d to 00 
in Lemma |8] part (a) is proved. Moveover, substituting c with 
and d with b^ in Lemma |S] part (b) is also proved. ■ 
These results can be easily generalized to the standard 
Nuttall Q— function. 

Corollary 3. Let /i, > and — v > 1 be fixed. Then the 
following assertions about the standard Nuttall Q-function, 
Qii.v (a, b), are true: 

(a) The function a i—s- {\/o., b) is log-concave on (0, 00) 
for b>0. 

(b) The function a ^ a'^/^K 2 fx(M-'^-^)/2) _ 
Qfj.,v iV^', b) is log-concave on (0, 00) for b > 0. 



VI. Bounds of the generalized Marcum and 
Nuttall Q-functions and their tightness 

A. Closed-form expressions of the generalized Marcum and 
Nuttall Q^functions with special order 

Recently, a closed-form expression of the generalized Mar- 
cum Q— function, Q^{a,b), was proposed for the case when 
V is an odd multiple of 0.5, given by [32. eq. (11)] 



Qy{a,b) 
-erfc I — !pr 

1 

1 



-erfc I — ^ 
2 V V2 

fe 



fe^" (-l)^(2g)! 



2q 

xE 

1=0 



1 



{ab) 



2q-ij\ 



a > 0,6 > 0. 



(38) 



C. The log-concavity of the Nuttall Q—function with respect 
to a 



X 



For the first step, we still consider the integral of 

/^2 (a;) on an interval (c, d) C (0, 00), like 



{m-v-1)/2 



for the generaUzed Marcum Q— function. The following result 
is true. 

d 



Lemma 8. The function a S ^ ^^^^fx^, ,,^ {^)dx is 



"^2(,y + l).a 

log-concave on (0,oo) for fixed /i — > 1, (c, d) C (0,oo) 
and v > 0. 



where erfc(-) is the complementary error function |49 eq. 
(7.1.2)]. For the case a — 0, the value of the generalized 
Marcum Q— function is 1321 eq. (12)] 



Q.(0,5) 



= erfc I —= 

.V2 



I'-l-S ^2fc+l 

E/ ofc-l $Z 



(-1)^ 



2- 



-(fc-9)!g!(2g + l) 



,(39) 



where v is an odd multiple of 0.5. 

More compact closed-form expressions for the case a > 
were derived in 1211 and 1331 . based on the recursion formula 
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of Q,(a, b) f? eq. (4.34)], given by 

= Q{b + a)+Q{b-a) 



n=l 



n-1 

E 

fc=0 



kl 



+ 0.5 e N, a > 0, (40) 

where (•). is the Pochhammer's symbol 1491 eq. (6.1.22)], 
Q{x) = ierfc (^-^^ and 

2 V V2 y 2 V V2 

, , y-1.5 / '=+0.5 



•\/27ra6 

Lr=0 



e 2 

fc=o 
(fc + r)! 
r\{k - ry.{2aby 



[(-l)'^e°'' + (-l)'=+ie 



fc+l„-ah] 



a>0,6>0. (41) 



We note that the main difference of ( I40> and ( I41> is that 
different formulas of the modified Bessel function of the first 
kind were used during the course of their derivations. For the 
case a — 0, the generalized Marcum Q— function can be also 
expressed as O eq. (4.71)] 



.(0,6) 



(42) 



Til.) ' 

where r(-, •) is the upper incomplete gamma function f49 eq. 
(6.5.3)]. 

If z/ is a positive integer, i42i reduces to |4. eq. (4.73)] 

Q.(0,5) = e-^gi^. (43) 
^-^ m! 

m— 

Moreover, we can derive a novel formula for Q,y{0, b), when 
ly is an odd multiple of 0.5. It is known that |49 eq. 8.356.2] 

r(a+ Lx) ar(a,x) +x"e""=, (44) 

and EH eq. 8.339.2] 

r{n + ^) = ^{2n-l)U, (45) 

where n is a positive integer and (2n+ 1)!! = 1 • 3 . . . {2n+ 1) 
1163 1 p. xliii]. Therefore, we can get the recursion formula of 
(5i/(0, b) for v odd multiple of 0.5 after some manipulations 



g,+i(0,6) = Q,(0,6) + e-- 
From the integral of i2Qi . we have 

b 



where 6 > 0. 



2 b' 



2iy 



TT (2;/)! 



Qo.5(0,6)=erfc( — ) = 2Q(6), 



(46) 



(47) 



Therefore, we obtain a new closed-form expression of 
Q^{0, b) for V odd multiple of 0.5 



Q^(0,5) =erfc 



V2 



e 2 



2 i^-l-S ^2/c+l 



E 



, (48) 



which is more compact than ( I39> . 

Similar result was proposed for the standard Nuttall 
Q— function in 1211 Theorem 1], given by 



(-l)"(2a)-"+^ 



E 

k=0 



(49) 



where a > 0, & > 0, /i > z^, m = /i+0.5 G N, n = i^+0.5 G N 
and the term „(a, b) is given by 



i=0 



+ (-1)' 



1— n— ^— It 



l + l {b + a^ 



[sgn(6 - a)] 7 ' 



(50) 



where 7(2;, x) = r(z)— r(z, x) is the lower incomplete gamma 
function |I49I eq. (6.5.2)], ( ' ) is the binomial coefficient |63 , p. 
xliii] and sgn(-) is the signum function. We note that the term 
Xj'jj „(a, b) can be further expressed by erfc(-) and exponential 
function with the help of (|^, (g^l and (gljl. 

The normalized Nuttall Q— function Q^.,y{a,b), for m = 
/I + 0.5 e N, n = J/ + 0.5 e N, can be evaluated by (21] 
Corollary 1] 



Qfi,i^[a,b) 
(-1)"2- 



T^E 



fc=0 



(51) 



where the term „(a, 6) is given by ( l50l l. 

B. Upper and lower bounds for the generalized Marcum and 
Nuttall Q— functions 

In 1 2 1 1 and 1 32 1 , upper and lower bounds for the generalized 
Marcum and Nuttall Q— functions were proposed by using the 
monotonicity of these functions. We find that the log-concavity 
of the generalized Marcum and Nuttall Q— functions can be 
used to establish even tighter bounds. 

1 ) Upper and lower bounds for the generalized Marcum 
Q~function with real order v: Although a lot of previous 
works only considered the generalized Marcum Q— function 
Q^{a,b) of integer order, the bounds for Qy(a,b) of non- 
integer order are desirable. Let [xj be the maximal integer 
less than or equal to x. Then, vi = \v + 0.5J + 0.5 is the 
minimal order that is larger than v and also an odd multiple 
of 0.5, V2 = Yv — O.5J + 0.5 is the maximal order that is less 
than or equal to v and is an odd multiple of 0.5. Since we have 
closed-form formula of Q,y{a, b) for ly that is an odd multiple 



IEEE TRANSACTIONS ON INFORMATION THEORY 



11 



of 0.5. The log-concavity of ly Qi/{a., b) on [1, oo), given in 
part (b) of Theorem |3] imphes one lower bound of Q^{a, b) 

Qu(.a,b) > Q,y-LBi{a,b) 

= Q.i(a,fo)"""^Q..(a,&)"^"", >^ > 1.5. (52) 

and also two more inequalities 

Qi^i{a,b) 

> g^(a,6)^^r^Q^,+i(a,6)^^^, 1, (53) 
Q^^ {a,b) 

> g,(a, 6) ^^Q,,„i(a, 6)^^, i/> 2.5. (54) 

After some simple algebraic manipulations, two upper 
bounds for Qu{a, b) are obtained, given by 

Qu{a,b) < Qu~UBi{a,b) 

= QuAa,bf''+^ IQu.+iia.by'-" , v>\, (55) 

and 

Qy(a,b) < Qu-UB2{a,b) 

= Q,, (a, fe)"""^+VQ.2-i(a, , v > 2-5. (56) 

Recall that in Theorem [2 we obtained that the function 
1-^ Qv{o-i b) is strictly increasing for v G (0, oo). Using this 
result, we can easily obtain that 

Qu2{a,b) < Qu^LBi{a,b) < Q^{a,b) 

< Qu-UBi{a-b) < Qy^{a,b), (57) 

where a > 0, 6 > and v > 1.5. This means that our bounds 
Qu-LBi{a,b) and Qu-UBi{a,b) are tighter than Q^^{a,b) 
and Qui{a,b) proposed in ||2T1 . ll32l . 

We can also use the log-concavity of the function ly 
Q^{a,b) — Q^{Q,b) on (0,oo), given in part (c) of Theorem 
|3] to get new bounds of Q^{a,b). We have the following 
inequalities 

QAa,b)-Q,{0,b) 

> [Q.^ia,b) - Q,jO,b)r-''' 

X [Q,, (a, b) - Q,, (0, 5)]"^"" , ly > 0.5, (58) 
QuAa,b)~Q,,{0,b) 

> [QAa,b)~Q,{0,b)]^T^ 

X [Q,,+i(a,6)-Q,,+i(0,5)]^r^, ly > 0, (59) 
Q,,(a,6)-Q,,(0,6) 

> [QAa,b)-Q,iO,b)]^^ 

X [Q,,_i(a,6)-Q,,_i(0,5)]^^, iy>1.5. (60) 

Therefore, we obtain another lower bound of {a,b), given 

by 

Qu{a,b) > Qy-LB2{a,b) 
= Q.(0, b) + [Q,, (a, b) - Q,, (0, b)]""-"' 

X [Q,, (a, b) - Q,, (0, b)p-'' , ly > 0.5, (61) 

and two upper bounds 



Qu{a,b) < Qy^uBz{a,b) 
= g.(0, b) + [Q,, {a, b) - Q,, (0, b)r-'+' 

I [g.i+i(a, b) ~ g,,+i(0, b)^'-'' , 1/ > 0, (62) 

and 

Q„{a,b) < Qu-UBi{a,b) 

= g.(o, b) + [g,, (a, b) - g,, (o, 6)]''-"=+' 

/[g,,_i(a,6)-g,,_i(0,6)]'^-''% z.>1.5. (63) 

We note that the log-concavity of the function v ^ 1 — 
Q^{a,b) can be also used to generate bounds of Q^{a,b), 
but the derived bounds are not reliable for large b. 

2) Upper and lower bounds for the Nuttall Q— function 
when — v > \ is an integer: Let /ii = [/i + 0.5j -I- 0.5 
is the minimal order that is larger than /i and also an odd 
multiple of 0.5, IJ.2 — La* ^ O-^J + 0.5 is the maximal order 
that is less than or equal to ji and is an odd multiple of 0.5. 
In Theorem |6] we obtained that the function v Qfi.u (oj b) 
is log-concave on [0,cx)) for ^ — v >1 fixed. Therefore, we 
can get a lower bound for the normalized Nuttall g— function 
when /X — > 1 is an integer, given by 

Q^L,u{a,b) > Qfj„u-LB{a,b) 
- Q^.,,uAa,br-''Q^,,,,{a,br-', > 0.5. (64) 

Moreover, we can get two upper bounds for Qf_i,,^{a,b) with 
/i — > 1 an integer after some simple manipulations, given 
by 

Qf^,i^ia,b) < Qf_,^„-^uBi{a,b) 
= Q^i,., (a, fo)'^^""+VQMi+i-i+iK ^>0, (65) 

and 

S^i,,.(a, 6) < Q^,u-UB2ia,b) 
= Sp2,.2(a,fe)'^^"^+VSp2-l,.2-l(«,^)''""^ i^>1.5. (66) 

In Theorem [2 we proved that the function i/ 1-^ 
Qp+,y4_i i,(a, 6) is strictly increasing. Therefore, we have that 

Q^2.f2(a, 6) < Qf^,i^-LB{a,b) < Qf_,,„{a,b) 

< Q^i,i^-UBi{a,b) < Qi_,^,^Aa,b), (67) 

where a,b > 0, > 0.5 and ji — v > 1 is an in- 
teger. This means the proposed bounds Q^.,y-LB{a,b) and 
Qij,,u-UBi{ci,b) are tighter than those given in |21|. 

These results can be also generalized to the standard Nuttall 
g— function, g^ j^(a, 6) with /i — > 1 an integer, with one 
lower bound given by 

g^. 1^(0, 6) > Qfi,iy-LB{a,b) 
= Q^„.Aa,br-'-Q^,^,,{a,br-'', iy>0.5, (68) 

and two upper bounds 

Qfj^,iy{a,b) < Qi_,,„-UBi{a,b) 
= Q^,.,Aa,bY'-''+' IQ^,+i,u,+i{a,br^-' , v>Q, (69) 
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Q, „„,(2.5.b)Q2_uB4(''''" 

Q, ,,.,(4,b) 




012345670123456 

b b 

(a) The bounds Q^„ifli(a, fe), Q„_j7fli(a, 6) and Qi,-uB2(a,b) (b) The bounds Qv-LB2(a,b), Qv-UB3{o.,b) and Q^-uB4,{a,h) 

Fig. 4. Numerical results for Qiy{a, b) and the proposed bounds versus b for a 6 {1, 2.5, 4} and u = 2. 





01 23456701 2345 

b 

(a) The bounds Q^_ifli(a, fe), Qy__j7fli(a, fe) and Qu-uB2(a,b) (b) The bounds Qv-LB2{a,b), Qv-UB3{a,b) and Q^-uBi{a.,b) 

Fig. 5. Numerical results for Qv{a, b) and the proposed bounds versus b for £ {2, 5, 8} and a = 2. 



and 

= Q^,,,,{a,br-'''+^/Q^,-i,,,-i{a,br-'\ iy>1.5. (70) 

We note that these bounds for the standard Nuttall Q— function 
Qpi,u{o-ib) are tighter than those given in |211 when a > 1 
by means of (I67> . Moreover, these bounds hold true when 
< a < 1, while the bounds given in 1211 do not. 

C. The tightness of the bounds 

1) The bounds of the generalized Marcut7i Q— function: For 
the bounds of the generalized Marcum Q— function Qi,{a, b), 
we first consider their tightness with respect to different values 
of the parameters. Fig. |4] shows the bounds of Qy{a,b) with 
different values of a, i.e. a E {1,2.5,4}, when v = 2; 
Fig. |5] shows the results for different values of i^, i.e. h' E 



{2,5,8}, when a — 2. We observe that the tightness of 
the bounds Qu-LBi{a,b), Qu-uBi{a,b) and Qu-UB2{a,b) 
improves as either a or v increases, and the tightness of 
the bounds Qu-LB2{a,b), Qu^UB?.{a,b) and Qu-UBi{a,b) 
improves as v increases, but worsens as a increases. There- 
fore, Qu^LBi{a,b), Qu~UBi{a,b) and Qv^uB2{a,b) are 
more suitable for large values of a, while Q,y^LB2{a,b), 
Qi^-UBsia, b) and Qv^uBiia, b) are proper for relative small 
a. However, since (5^(0,6) becomes very small as b grows, 
these two groups of bounds tends to be equal for large b. 

The numerical results of our bounds Qi^-LB2{o,,b), 
Qv-UB7i{a,b) and Q^-ijB4{o,,b) with non-integer order i/ 
are shown in Fig. |6l We can see that our new bounds are 
tighter than the bounds proposed in Uf I . It is interesting 
to see that the upper bound Qy-uBz{o-ib) is tighter than 
Qu-UBAio-^b) when v = 5.1, but Qv^uBiio-^b) can be 
tighter than Qu-UB7i{o-,b) when v — 1.8. This is because 
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vi — v — 0.4 < v—1'2 ^ 0.6 when v — 5.1, and the inequality 
in \59\ can be tighter than that in \6Q\ . When v ~ 1.8, 
vi — V = 0.7 > V — V2 — 0.3 and the inequality in ( I59> 
is looser than that in (I60> . Similar results can be also found 
for the bounds Qu-UBi{o-ib) and Qi/-uB2{a,b). Therefore, 
we can choose the tighter upper and lower bounds proposed 
in this paper according to the value of a and the decimal value 
of ly. 

Next, we compare the proposed bounds with other existing 
bounds when u is integer Since most of the existing bounds 
for the generalized Marcum Q— function of integer v order are 
valid for only either 6 > a or 6 < a, we choose to show the 
comparisons case by case. 

For the case b > a, the existing lower bounds include LBl- 
AT in [35,, the first line in eq. (18)]', Q„i-o.5(a,fc) in iSS 
eq. (11) and (14)], GLBml-KL in OH eq. (6)], LBl-BS in 
FT eq. (4) and (8)] and LBl-B in f4T^ eq. (8) and (15)]. The 
existing upper bounds include UBl-SA in [34. eq. (8)], UBl- 
AT in 135 eq.(17)], Q„+o.5(a,6) in M eq. (11) and (14)], 
GUBml-KL in f3F eq. (5)], UBl-BS in (4? eq. (5) and (9)] 
and UBl-B in |42 eq. (9) and (16)]. Fig.|7lshows the results 
for the case 6 > a = 1.5 and = 2 in a logarithmic scale. We 
choose relative small values of v and a in order to facilitate 
the recognition of our bounds from the exact value. Even in 
this case, our new bounds are shown to be much tighter than 
the other bounds in the literature. For larger values of v and 
a, our numerical results show that our bounds are very tight. 

For the case b < a, the existing lower bounds include LB2- 
SA in 1 34 eq. (12)], LBl-AT in |35 the first fine in eq. 
(18)], LB2-AT in |35 eq. (20)], LB3-AT in |33 eq. (21)], 
Q^_o.5(a,6) in |32, eq. (11) and (14)], GLBm2-KL in EH 
eq. (9)], LB2-BS in |41 eq. (11) and (17)] and LB2-B in 
1421 eq. (17) and (23)]. The existing upper bounds include 
Qm+o.bia^b) in 02 eq. (11) and (14)], GUBm2-KL in 01 
eq. (8)], UB2-BS in fTT" eq. (11) and (17)] and UB2-B in 
1421 eq. (18) and (24)]. The numerical results for the case 
b < a = 4 and ly = 4 are illuminated in Fig. |8l Our bounds 
Qi^-LBi{o.,b) and Qu-uBi{a,b) are much tighter than the 
other bounds in most of the cases, but the Qi,~uB2{a, b) can 
be looser than the existing bounds. 

Let the absolute relative error of a bound be e% ~ 100% x 
|bound-Q (a,b)| ^ If J/ = 4 and fl = 4, the maximal absolute 

relative errors of our bounds Qy-L_Bi(a, 6), Qiy-uBi{a,b) 
and Qu-UB2{a,b) are 0.4437%, 1.3077% and 1.3799%, re- 
spectively, for all range of b. If v = Q and a = 6, the 
the maximal absolute relative errors of our bounds reduce to 
0.2157%, 0.6420% and 0.6581%, respectively. Our numerical 
results suggest that the absolute relative errors of the proposed 
bounds are less than 5% in most of the cases. 

Moreover, the relative errors of all our bounds converge to 
as 6 — > cx). This can be explained simply by the asymptotic 
formula of the generalized Marcum Q-function when b 
oo, given in @. From @, we can find that Q^{a, b) inclines 
to log-linear for very large b and the inequalities i52i . i53i 
and ( I54t tend to be equal. Hence, the bounds Qi,-LBi{a,b), 

'There is a mistake in tlie formula given in I35il . For a correct version, the 
readers are refen'ed to equation 1431 . 
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b 

Fig. 6. Numerical results for the bounds of Q^(a, b) with non-integer order, 
where u = 1.8,5.1 and a = 1.8. 'x': exact. Dashed line: previous bounds. 
Solid line: some of our new bounds including Qv-LB2 (ij b), Qu—uB:i{o., b) 
and Q^_(7S4(a, b). 




2 2.5 3 3.5 

b 



Fig. 7. Numerical results for Q^(a, b) and its upper and lower bounds versus 
b for the case b > a = 1.5 and u = 2. 'x': exact. Dashed line: previous 
bounds. Solid line: some of our new bounds, including Q^— Lsi(a, fc), 
Qiy-UBlia,b) and Q„-uB2{a,b). 



Qu-UBi{o-;b) and Qu-UB2{o-,b) converge to the exact value 
as & ^ cx). Since (5,^(0,6) — > as 6 — > oo, we obtain that 
the bounds Qu-LB2{a,b), Q,y-uB3{a,b) and Qu-UBi{a,b) 
possess the same property. 

On the other hand, it was proved that the relative errors 
of the bounds Qi/-o.5(aj b) and Qu+o.bia-, b) does not tend to 
zero as b approaches infinity [41 . And the bounds proposed 
in 1411 . whose relative errors tend to zero when 6 ^ oo, hold 
true for only 6 > a. To the extent of the authors' knowledge, 
our bounds are the first bounds with such tightness on the 
whole region of 6 > 0, even in terms of relative errors. 

2} The bounds of the Nuttall Q~function: We compare 
the proposed bounds of the normalized and standard Nuttall 
Q— function of the order /x, i/ > with other existing bounds. 

The existing lower and upper bounds for the normalized 
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Fig. 8. Numerical results for Q,y(a, b) and its upper and lower bounds versus 
b for the case 6 < a = 4 and u = 4. 'x': exact. Dashed line: previous bounds. 
Solid line: some of our new bounds including Qv-LBi (Oj b), Qu—UB\ ifl, b) 
and Q„-uB2(a., b). 
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b 

Fig. 9. Numerical results for the normalized Nuttall Q— function Sfi,i/(a, b) 
and its bounds versus b for different values of a, i.e. a = 1, 3, when ^ = 4 
and V = 2. 'x': exact. Dashed line: previous bounds. Solid line: our new 
bounds including Q^,„-LB{a,b), Q^,v-U Bi{a.,b) and Qfj.,v-u B2{a,b). 



Nuttall Q-function Q^,i,(a, 6) are Ql^Jq.s.Mo s ^) 
2mo.5,Mo.5("'^) ED ^1- 19], respectively, when ^,>v + l 
and V > 1. Fig. |9] shows the numerical results for different 
values of a, i.e. a e {1,3}, when = 4 and = 2; Fig. ^| 
shows the numerical results for different values of v and ^, 
i.e. V G {1, 4} and ji = v + 2, when a = 2. We find that the 
tightness of our bounds Qf_i,u-LB{a,b), Qf_L,u-UBi{a,b) and 
Qfj..iy-UB2{a,b) improves as either a increases or i/ increases 
with fi — ly fixed. This result is expected for Q^^i,{a, b), since 
it holds true for the special case of the generalized Marcum 
Q— function Qu{a, b). As we have proved in (I67> . the proposed 
bounds Qf_i,v-LB{a,b) and Qf_i^v-uBi{a,b) are tighter than 
the bounds given in 1211 . 

The numerical results for the case that v is not an integer 
and ^ — V > 1 is an integer are shown in Fig. ^2 We 




1 1.5 2 2.5 3 3.5 4 



b 

Fig. 10. Numerical results for the normalized Nuttall Q— function Q^,i^(a, b) 
and its bounds versus b for different values of u and /i with fj, — u fixed, 
i.e. 1/ = 1,4 and fi — i/ = 2, when a = 2. 'x': exact. Dashed line: 
previous bounds. Solid line: our new bounds including Q^.i/— ls (fi, fe), 
Sji,i/-(7fli(a, f)) and Q^^i,-uB2{a, b). 




Fig. 11. Numerical results for the bounds of the normalized Nuttall 
Q— function Qf_i,u(o.,b) with non-integer order, where fi — u = 2, u = 
1.7,5.2 and a = 1. 'x': exact. Dashed line: previous bounds. Solid 
line: our new bounds including Q^,^_lb (a, b), Q^i, ly-c/si (a, fe) and 
Q/i,i/-(7S2(a, b). 

can see that the tighter ones of our new bounds are tighter 
than the bounds proposed in |21| when the orders /i and 
1/ are not integer Moreover, we find that the upper bound 
Q^i,v~UBi{a,b) is tighter than Q ^^^u-UB2{a,b) when v = 5.2, 
while Qf_,,„-uB2{a,b) can be tighter than Qp,,y-;7si(a, 6) 
when v — 1.7. This observation is quite similar with the 
case of generalized Marcum Q— function. It means that we 
can choose the tighter upper bounds of the normalized Nuttall 
Q— function according to the decimal value of v. 

The numerical results for the exact value and bounds of 
standard Nuttall Q— function are shown in Fig. [21 When a > 
1, the proposed bounds are tighter than the bounds given in 
|21 1. When a < 1, our proposed bounds hold true, while the 
bounds of 1211 do not. The bound Q^.i/-uBi{a,b) is also 
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b 

Fig. 12. Numerical results for the standai'd Nuttall Q— function Qfj,,v{a, b) 
and its bounds versus b for for different values of a, i.e. a = 0.5,1.5, 
when u = 3, fi = I. 'x': exact. Dashed line: previous bounds. Solid line: 
some of our new bounds including Q^.^— lb (a, fe), Qn^^—uBli'^yb) and 
Qij,,,^-UB2ia, b). 

tighter than Qf^^u-UB2{a, b), and we can use Qfj.^u-UBi{a, b) 
in the appHcations. The other properties for the bounds of 
normalized Nuttall Q— function maintain for the bounds of 
standard Nuttall Q— function in terms of relative errors. 

Now we consider the tightness of our bounds for the Nuttall 
Q— function on the whole region of b E (0, oo). For a bit 
larger parameters, when /i = 7, = 4 and a = 4, the 
maximal absolute relative errors of our bounds Q ^,v-lb{o-i b), 
Q^,M-UBl{a,b) and Q^,,^^uB2{a,b) are 0.4134%, 1.2216% 
and 1.2802%, respectively, for all range of b. When fi — 9, 
V — Q and a = 6, the maximal absolute relative errors 
of our bounds reduce to 0.2079%, 0.6190% and 0.6334%, 
respectively. Our numerical results suggest that the absolute 
relative errors of the proposed bounds are less than 5% in 
most of the cases. 

Using the same method given in |4 p. 81], we can get the 
asymptotic formula of the normalized Nuttall Q-function 

Q^^^{a,b)^-^Q{b-a). (71) 
For the standard Nuttall Q-function, we have 

Q^^,{a,b)^ —^Q{b-a). (72) 

It implies that the relative errors of our bounds of the nor- 
malized and standard Nuttall Q— functions also converge to 
as 6 ^ oo. Moreover, using M\\ it is quite simple to 
show that the relative errors of the bounds Qy^\a r,,Li^Jo sC*^' ^) 
and Q[^t]o 5,[iy]o 5 1*^' ^) ^'^^ ^^^'^ '■° ^ approaches 

infinity. To the extent of the authors' knowledge, our bounds 
for the normalized and standard Nuttall Q— functions are the 
first bounds with such tightness in terms of relative errors on 
the whole region of 6 > 0. 

VII. Applications 

The proposed bounds for the generalized Marcum and 
Nuttall Q— functions have been shown to be quite tight in 



the previous section. They involve only exponential function 
and the erfc function, and therefore can be computed very 
efficiently. These results can be applied to the performance 
analysis of various wireless communication systems operat- 
ing over fading channels. Some application examples of our 
proposed bounds are given as follows: 

One interesting example occurs when bounding the outage 
probability of wireless communication systems which are 
both interference- and power-limited, where the Rayleigh/ 
Nakagami faded desired signals are subject to independent 
and identically distributed (i.i.d.) Rician faded interferers |'64J 
eq. (18) and (26)]. Another similar application is to bound 
the outage probability of maximal ratio combining (MRC) 
in the presence of independent but not necessarily identically 
distributed Rayleigh faded co-channel interference, when the 
received signal at every antenna experiences i.i.d. Ricean 
fading flT eq. (21), (27) and (28)]. The techniques in these 
two examples apply directly in the outage probability analysis 
of cognitive radio (CR) system and wireless sensor networks 
(WSN) f67l. 

Our proposed results can also be applied to evaluate the 
average error probability of digital communication systems 
operating over slow-fading channels. In |25l, the authors ap- 
proximated the average error probability of MRC multichannel 
reception by using a piecewise polynomial approximation 
method. Their results require to compute the partial (truncated) 
general moment of generalized Rayleigh or Ricean distribu- 
tion, which can be simply represented by generalized Marcum 
and Nuttall functions. Actually, comparing with the infinite 
series formulation in |25 eq. (17)] for the case s > 0, one 
would prefer to use the result 

^) (1^),, 73, 

where is subject to a distribution of cr^x^s^/o-^- '^^'■^ 
that the order of generalized Rayleigh or Ricean distribution 
can be not integer Therefore, the evaluation of non-integer 
order generalized Marcum and Nuttall Q— functions is needed. 

Finally, our results can be used to extract the log-likelihood 
ratio for the decoding of turbo or low-density parity check 
(LDPC) codes for differential phase-shift keying (DPSK) 
signals L65J . 

VIII. Conclusion 

In this paper, we provided a comprehensive study of the 
monotonicity and log-concavity properties for the generalized 
Marcum and Nuttall Q— functions. Tight upper and lower 
bounds for the generalized Marcum and Nutall Q— functions 
have been obtained by using the log-concavity of these 
functions. If the bounds are chosen based on the values of 
the parameters, the proposed bounds are tighter than the 
existing bounds in the literature in most of the cases. Our 
proposed bounds are tight in terms of relative errors for all 
range of b. We have proved that the relative errors of our 
proposed bounds converge to as cx). The numerical 

results show that the absolute relative errors of the proposed 
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bounds are less than 5% in most of the cases. To the extent 
of the authors' knowledge, our bounds for the generalized 
Marcum and Nuttall Q— functions are the first bounds with 
such tightness in terms of relative errors on the whole region 
of 6 G (0,cx)). Some applications of the proposed theoretical 
results have been also provided. Further research directions 
include to prove the conjectures that were provided in Section 

EYl 
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Appendix A 
Proof of Lemma[T] 

It is known that Ex{EYiY\X = x)) = EyiY) and P{Y e 
S) — Ey{Cs), where £s is the indicator function of the set 
S. Therefore, 



P(r eS) = Ex{P{Y e s\x = x)). 



(74) 



Since X and Y are non-negative independent random vari- 
ables, we have 

Fx+Y{t) = P{X + Y <t) 
^Ex{P{x + Y <t)\X ^x)) 



Frit - x)fx{,x)dx. 



(75) 



The differential of Fx+Y{t) is 
fx+Y{t)dt 

= FY{Q)fx{t)dt + I dEvit - x)fxix)dx, (76) 
Jo 

where we have used the definition of Riemann-Stieltjes inte- 
grals. Since Ey{0) > 0, we have 

fx+Y{t)dt> f dEY{t-x)fx{x)dx. ill) 

Jo 

This in turn implies that 

poo 

/ 9{t)fx+Y(t)dt 

Jb 

> 9{t) ^ dEY{t ~ x)fx{x)dx 
= / / g{t)dFY{t-x)] fx{x)dx 

Ja \Jmax{x,b} j 

(/ 9{t)dFY{t^x)\j]x(x)dx 

POO / POO 



fO \Jb-x 

POO / POO 



g{t)dEY{t-x)j fx{x)dx 
g{x + y)dFY{y)^ fx{x)dx 

g{x + y)dEY{y)) fx{x)dx. (78) 



Here, we note that the first term and the second term in above 
result are exactly the integrations on region A and B (see Fig. 
[0, respectively. Since Fy(0) < 1, we know dEyiy) > on 
a subset of (0,oo) with non-zero measure. This indicates the 
following inequality 

g{x + y)dEY{y) ] fx{x)dx > 0. (79) 

/O \Jb-x J 

On the other hand, we know that 

roo / roc \ 

i^J^ g{x + y)dEY{y)j fx{x)dx 

/OO / p oo \ 
^ dFY{y)^g{x)!x[x)dx 



g{x)fx{x)dx. 



With this, the proof is complete. 



(80) 



Appendix B 
Proof of Theorem^] 

(a) The proposed result is equivalent with the log-concavity 
of Q^{0,Vb). We have 



f^2jx)dx 



fxl {^)^[b,oo){x)dx, 



(81) 



where Cip,oo) is the indicator function for the interval [6, oo). 
We can easily prove that ^ C\p ,^){x-\-y) is SRR2 for 

X, y > and 6 > from the definition. Using the reproductive 
property of the central chi-square distribution one has 



/v2 ix) — /v2 ix) * f^2 (x) 



(82) 



Applying parts (a) and (b) of Lemma |5] for the pdf of 
(central) chi-square distribution given in ilQ\ we get that 
(x, ly) 1-^ f^2 (x) is TP2 for x,i^ > 0. Then, from Lemma 0] 
we can easily obtain that (^^^^^^(O, Vb) is SRR2 for vi > 
and 1^2 > 0. This implies that ly 1-^ (3^(0, Vb) is log-concave 
for > due to part (d) of Lemma |5] Substituting b with 6^, 
part (a) is proved. 

(b) In view of (|8} and (I74> . we get that 

Q^+i (^^/a, Vbj 

= P{xlu+2,a>b) 

= P{xlu + xla>b) 

^E^2(P{xl,+x>b\xl,^x)) 



.fxlM)P(X2,a >b-x)dx, 



(83) 



where 



Pixla>b-X) 



Qi {y/a, \Jb — x) , < x < &, 
1, x>b. 



(84) 
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We first note that x ^ PiX2 a — ^) continuous on 
(— cx3,cx3). Recently, we proved \44] that b i-^ Qy{a,\/b) is 
strictly log-concave on (0, oo) when a > and v > \. Using 
the same proof process as in |44|, we can also obtain that 
h ^ Quia, Vb) is log-concave on (0, oo) for = 1. In view 
of part (e) of Lemma |5j we have that x i— > Qi (a, \/b — x) 
is log-concave on (— oo, b) for a > 0. Therefore, it is 
easy to prove that x i-^ P{X2 a — b — x) is log-concave 
on (—00,00), since x t-^ Qi{a, \Jb — x) is increasing and 
x ^ PiX2 a ^ b — x) is continuous on (—00,00). Hence, 
{x,y) ^ P{x2,a >b~{x + y)) is SRR2 due to pait (d) of 
Lemma |5] In view of ( I82> and the fact that {x, v) ^ f^2 (x) 
is TP2 for x,i^ > 0, we can use Lemma@]to (I83> and obtain 
that {1^1,1^2) ^ Qyi+u2+i {-/a, \/bj is SRR2 for i/i, 1/2 > 0. 

This implies that ly Q^, y\/a^ Vb^ is log-concave for > 1 
and a,b > due to part (d) of Lemma|5] Since 1^ 1-^ f^2 ^ is 

continuous for ly > 0, the function v i—i- (y^a, \fbj is also 
continuous at i/ = 1 for a,b > 0. Hence, we can prove that 

1^ Qu (^/o-i is log-concave for v > 1. Substituting a 
with and b with b^, the asserted result is proved. 

(c) The function Qui^/a, Vb) can be also represented by 
using the non-central chi-square distribution with zero degrees 
of freedom, as follows 



We have 



Qu i^V^, Vbj 

'P{xl,a>b) 

-Pixl.+Xla>b) 

:E^2^ {P{x + xla>b\xl.^^)) 

i 

fxlA^^^(^o,a > b-x)dx. 



(85) 



All the same, x P{Xx) a — ^) has a discontinuous point 
at a; = 0, so it is not log-concave on the whole (—00,00). 
However, by using (I15> we get 



Qu ( \fa, Vb 



fxlM)P(X^,a > b-x)dx 



fxL (^)^^ + / 4i„ (2^)-P(xL >b- x)dx 



= Qu (0, Vb)+ f f^2^ ix)PixL > b - x)dx, (86) 



where 



, , , /"^ 1 /a\i/2 , , 

P {Xla ^ ^) = / 2 \x) "^"^ ^^^^ 

for all 6 > 0. Now, let us consider the normalized modified 
Bessel function of the first kind -fu : (0, 00) 1-^ (1, 00), defined 
by 



P {xla > b) 
00 -j^ 

—^i{ax)ae~ 
b 4 
00 ^ 

-'ji{ax)ae' 



'2° da; 



-c 



[6,00) 



{x)dx 



^ji{ax)ae ^ C[o.oo){x - b)dx. (89) 

It is known that the function x t-^ lu{x) is log-concave if 
V > ~1 1661 Theorem 2.2]. So x ^ ^i{ax) is also log- 
concave for a > 0. On the other hand, since the set [0, 00) is 
convex, the indicator function x ^^[0,00) (2;) is log-concave. 
Hence, the function (a;, b) ^ £[o,oo)(2^~^) is also log-concave 
due to part (e) of Lemma |5l Therefore, the integrand of ( I89> is 
log-concave in {x,b), as a product of log-concave functions. 

It is known that if (a;,?;) 1-^ / is log-concave, then 
g{x) = J f{x, y)dy is a log-concave function of x 1471 p. 106]. 
Applying this to ( I89> . we obtain that b ^ P [xo a ^ ^) is log- 
concave on (0, 00). It follows that the function x 1^ P{Xs> a ^ 
b — x) is also log-concave on (0, 6), due to part (e) of Lemma 
|5] Then, by the same method used in part (a) and (b), we can 
prove that the function 1/ 1— > /p (x)P(xo_a ^b— x)dx is 
log-concave on (0,oo). Hence, Qu{\/a,Vb)~-Qu{0,Vb) 
is log-concave on (0,oo). Substituting a with and b with 
b^, the asserted result is proved. 

(d) Similar to J85> . we obtain that 

1-Qu (0,Vb 



(90) 



Pixl<b) 



where £(o.b) is the indicator function for the interval (0, b). 
We can easily prove that (a;, y) 1-^ £(q ;,)(x + y) is SRR2 for 
x,y > and 6 > from the definition. Thus, using the same 
method as in the proof of part (a) and (b), we can prove that 
the function ly 1-^ 1 — Qu (0, b) is log-concave on (0, 00) for 
each 6 > 0. 

(e) Similar to ( I85t . we obtain that 



1 - Qu+i yVa, Vbj 
Pixlu+2.a<b) 

P{xlu+Xl.a<b) 

-■E^2jP{xla+X<b\xlu^^)) 
b 



= / f^2jx)Pixia<b-x)da 



f^2^ (x) (l - Qi (^^/a, Vb - x^^ da 



7,(x) = 2^r(z/ + l)x-'^/2j^(Vi) 



(88) 



dx. (91) 



We have proved the function b t-^ 1 — Qi(a, Vb) is log- 
concave on (0, 00) for all a > in Theorem|2] It follows that 
the function x t-^ 1 — Qi(a, Vb — x) is log-concave on (0, b), 
due to part (e) of Lemma |5] Thus, using the same method as 
in the proof of part (a) and (b), we can prove that the function 
ly 1-^ 1 — Qi, (a, b) is log-concave on [1, 00) for each a > 0, 
& > 0, and with this the proof of this theorem is complete. 
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We know that 



Appendix C 
Proof of Lemma|^ 



f^2 (x)dx 

= P{xlu.a^+a, e (C,d)) 

= P {Xlu,a, + Xla, e (C, d)) 

^ f P{x + xL. e (c, d)) „ (92) 

where we used that ^ (a; + Xo 02 ^ ("^i '^)) = ^ for all x > d. 
Recall that |46 p. 437] 

^(^ + X0,a2 € {C,d)) 



Appendix D 
Proof of Theorem|5] 

(a) We still use the normalized modified Bessel function of 
the first kind 7,^, which is defined in ( I88> , to represent the 
Nuttall Q— function. When a > 0, we have that 



00 ^(p-l)/2 



2a''/2 



-{x+a}/2 



I, J (Vox) dx 



(M+f-l)/2 



a; 



(m+i'-1)/2 

^''^''^^2-+ir(i. + i)' 



E 



-"^/2p (a; + xl e (c, d)) , (93) The function x 



for all Q < X < d. Moreover, for all j >Q integer P{x + X2j •= 
(c, d)) can be further expressed as 

/xij (^)'C{i;+te(c,d),x,t>o}(a;, (94) 

where C{xj^t£(c,d),x,t>o}{xi't) is the indicator function of the 
set S = {(x, t)\x + t ^ (c, d),x,t> 0}. 

Observe that (x, t) ^ C{^+te{c,d),x,t>o} is 5'i?i?2 for x, t > 
from the definition. Recall that from the proof of part (a) 
of Theorem |3] we already know that {t,j) 1-^ fx?, i^) i^ 
for t > and j > integer Applying Lemma |3]^or ( I94L we 
obtain that (x, j) P{x + X2j ^ (^i c')) is SRR2- From parts 
(a) and (b) of Lemma |5] we have that 

0-,a2)--^^e"-/2 (95) 



(m+i/-1)/2 

>M 2.+ir(^ + ^) e-(--+°)/^/:[o,oo)(a: - 6)dx.(97) 

Since the function x 1— > 7,^(2;) is log-concave if > — 1 
l66l Theorem 2.2], x ^ 7i/(ax) is log-concave for a > 0. 

2''+^r(v+i) is also log-concave 

when /i + > 1. On the other hand, since the set [0, 00) is 
convex, the indicator function x ^ C.{n.x!){x) is log-concave. 
Hence, the function (x, 6) ^ C[o.oo){x~'b) is also log-concave 
due to part (e) of Lemma |5] Therefore, the integrand of ( I97> 
is log-concave in (x, b) as the product of several log-concave 
functions when /i + > 1 and a > 0. Recall that if {x,y) 1— > 
/ is log-concave, then g{x) — J f{x,y)dy is a log-concave 
function of x |47 p. 106]. Applying this to (|^} and changing 
a with a^, part (a) is proved. 

(b) Since b t-^ Qp,j/(a, is log-concave and decreasing, 
the function g{b) = b^ is convex, then the composite function 
6 (/ o g){b) — Q^^u {a,b) is log-concave too |47 p. 84], 
and part (b) is proved. 

(c) When a > 0, i;^,,^^,,,^!^^— i) - Q^,.(^A,\/6) 
can be expressed as 



is TP2 for all j > integer. Using again Lemma |3] for J93> . 
we obtain that (02, x) 1-^ P(x + Xo,a2 ^ i^' d)) is SRR2 for 
02 > and < X < d. It is known that (x, ai) 1-^ f^2 (.x) 
is TPoo on X and ai with constant 1/ > fST] p. 110]. Then 
by another application of Lemma |3] for ( I92> . we deduce that 



^2'^+ir(i/ + l 



2:+ir(^ + 1) ^'^^+"^^'^(-00,.] {x)dx 

{p+u-l)/2 



(ai, 02) 



(x)dx 



(96) 



is SRR2 for ai and 02 > 0. This implies that indeed the 

d 

function a ^ j /^p (x)dx is log-concave for a > due to 

c _ 

part (d) of Lemma |5] Since the pdf of non-central chi-square 
distribution given in (0 is replaced by the pdf of central chi- 
square distribution given in ( I10> when a ^ 0, the function 

d 

a I— > / (x)dx is continuous at a = 0. Hence it is also is 

c 

log-concave for a > 0, which completes the proof. 



Therefore, using the same method as in the proof of part (a), 
we can prove easily the asserted result, 
(d) From (|2}, we obtain that 



-^X2(,. + l),a(^ 



fl—jy— 1 



a" Jo 



^I^{at)dt 



— / t^'t 



■I,{at)q_^^o]{x-b)dt. (99) 



We have proved that function x 1— > xl^{x) is log-concave 
for each v > 1/2 1441 Proposition 2.1]. Since {x,b) 
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'C(-oo.o](2; — b) is log-concave as proved in part (a), the 
integrand in ( I99l l is also log-concave in x and b for all /i > 1, 
ly > 1/2. Hence, using the same method in the proof of part 
(a), the asserted result follows. 

Appendix E 
Proof of Lemma0 

We know that 



^(M-.-i)/2 ( r y^^^ j^^^ _ i)^^ j 



\ •/ maa;{c,i} 







/ 7naa;{c— i,0} ' / 



Since / 2 (x) = for a; < 0, we have 



X2„ 



x(/ {x + tY''--'-'^^^f^2 {x)C[,^d){x + t)dx)dt.(m) 



When ^—v is fixed and no smaller than 1, a; 1-^ 2;(a'-i'-i)/2 jg 
log-concave. We also have that the function x 1— > C^c.d){x) is 
log-concave, since [c, d) is a convex set. Then, from part (e) of 
Lemma|5] we can conclude that {x,t) i— > (x + t)'^^^'^^^^/^ and 
(x, 1-^ C\c^d){x + t) are both log-concave. We know that the 
function a; ^ f^2 (x) is log-concave on (0, 00) for all a > 

|^|. Therefore, {x,t) ^ (x+t)^''-''-^^/^ f^2 ^{x)Cic,d){x+t) 

is log-concave as a product of log-concave functions. Hence, 

00 

we have that 5 (t) = / (x + i)*''"'""^)/^^! {x)£[c,d){x+t)dx 

is log-concave in t, and thus (ii,t2) '-^ 9{ti + ^2) is SRR2 
due to part (d) of Lemma |5] 

On the other hand it is known that 



f^2 (x) — f^2 (x) * f^2 (x) 



(102) 



and (x, J^) I— > ^2 (x) is for x,iy > 0. Then, 

applying Lemma in (I101> . we obtain that {1/1,^2) ^ 



i: 



2,(M-i^l-i^2-l)/2/- (2,^3, js 5-^^ foj. > 

fd^(M-l^-l)/2. is 

Jc ■'X2(„ + l),a^ ' 

> 1, due to part 



0, which implies that u 
log-concave for > and fixed fi 
(d) of Lemma [S] Moreover, it can be easily verified that 



u I— > r"x'^ ^)/-^f^2 (x)dx is continuous at = 0. 

Therefore, ly 1^ f x*^^^''^^-'/^ f^2 (x)dx is log-concave 
on [0, cx)) for fixed /i — > 1, and the proof is complete. 



Appendix F 
Proof of Lemma[S] 

Similarly with the argument of Lemma we can get 



X 



^2(v + l),ai+a2 



(x)dx 



it) 



f: 



X2 



(x)(x + tY^'-'^-^^/^Cicdjix + t)dx^ dt. (103) 

When fM — v is fixed and no smaller than 1, x i-^ 2;(a'-''-i)/2 
is log-concave. We also have that the function x 1-^ C[c,d){x) 
is log-concave, since [c, d) is a convex set. Then, the function 
X x^'^~'^~^'>^^C[c.d)ix) is also log-concave as a product of 
two log-concave functions. From part (d) of Lemma|5] we can 
conclude that the function (x,t) i-> (x+i)(^"''"i'/2£[c,rf)(x+ 
t) is SRR2 for X, t > 0. It is known that (x, i— > f^2 (x) 
is TPoo on X and 02 |E1 p. 110]. Using Lemma |3] we obtain 
that the function 

/■oo 

9{t,a2)= 42 {x){x + t)^^—^'>/^q,,a){x + t)dx 
Jo '"'^ 

(104) 

is SRR2 for 02 > 0, and 



^(m-!^-1)/2 



^2(,/ + l),ai+a2 



(x)(ix 



/2 (t)g{t,a2)dt. 



(105) 



We know that (i,ai) i-^ /^2 (t) is also TPqo 
for fli > when ly > 0. Therefore, by an- 
other application of Lemma |3l in (I105> . we get that 
(01,02) ^ f''x(^-''-i)/2f 2 {x)dx is 5'i?i?2 for 

z^, ai,a2 > 0. This implies that the function a > 
x^f^^"^^^^^ f^2 ^ {x)dx is log-concave on (0, 00) for 
> and /i — > 1 fixed, due to part (d) of Lemma |5] Since 
v 1^ f x^f^^"^^')^^ f^2 (x)dx is continuous at = 0, 
we can prove that the log-concavity holds true for the case 
ly = Q. Thus, the proof is complete. 
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